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ABSTRACT 

We present observations and analysis of an unusual [CII] emission line in the very lumi¬ 
nous QSO SDSS J155426.16+193703.0 at z ~ 4.6. The line is extremely broad (FWHM 
735 km s -1 ) and seems to have a flat-topped or double-peaked line profile. A velocity map of 
the line shows a gradient across the source that indicates large-scale rotation of star-forming 
gas. Together, the velocity map and line profile suggest the presence of a massive rotating disc 
with a dynamical mass Md yn > 5 x 10 10 M@. Using the assumption of a rotating disc ori¬ 
gin, we employ an empirical relation between galaxy disc circular velocity and bulge velocity 
dispersion (cr) to estimate that a > 310 km s _1 , subject to a correction for the unknown disc 
inclination. This result implies that this source is consistent with the local M-a relation, or 
offset at most by an order of magnitude in black hole mass. In contrast, the assumption of 
a bulge origin for the [CII] emission line would lead to a conclusion that the black hole is 
nearly two orders of magnitude more massive than predicted by the M-a relation, similar to 
previous findings for other high-redshift QSOs. As disc rotation may be a common origin for 
[CII] emission at high redshifts, these results stress that careful consideration of dynamical 
origins is required when using observations of this line to derive properties of high-redshift 
galaxies. 

Keywords: quasars: individual: SDSS J155426.16+193703.0-quasars: supermassive black 
holes - quasars: emission lines - galaxies: kinematics and dynamics - radio lines: galaxies 


1 INTRODUCTION 

With the recent inauguration of the Atacama Large Millimeter/Sub¬ 
millimeter Array (ALMA) telescope, the astronomical litera¬ 
ture has shown a surge in studies of the 157.74 pm forbidden 
2 ? 3/2 2 Pi /2 fine structure line ([CII]) of ionized carbon. 

This line arises from photo-dissociation regions in the cold neu¬ 
tral interstellar medium (ISM), marking it as a tracer of cool gas 
and star formation (for a recent review, see lCarilli & Walteiil2013l) . 
Furthermore, [CII] is typically the brightest emission line in the 
ISM: in ultra-luminous infrared galaxies (ULIRGs), this single line 
can account for ~ 0.1- 1 per cent of the entire far-infrared (FIR) 
luminosity l iLuhman et all19981) . Therefore, [CII] serves as an in¬ 
credibly powerful tool for probing the dynamics of star formation 
(SF) activity in galaxies. This statement is especially true at high 
redshifts owing to the high luminosity of [CII] compared to other 
galaxy emission lines, as well as to the high angular-resolution of 
ALMA. 

The [CII] line will be particularly valuable for observing 
quasi-stellar objects (QSOs, also commonly referred to as quasars), 
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because it is one of few tools to directly observe the host galaxies 
of these luminous sources. QSOs are the signposts of actively ac¬ 
creting super-massive black holes (SMBHs), and at most frequen¬ 
cies they vastly outshine their host galaxies. The SMBH’s accretion 
disc radiates strongly in the optical and ultraviolet (UV) regimes; 
the disc’s corona up-scatters disc photons into the X-ray regime; 
and hot dust close to the SMBH radiates thermally with a peak in 
the mid-infrared (MIR). However, in the FIR and sub-millimetre 
regimes, cool dust heated by star formation processes in the host 
galaxy has a substantial contribution to the emission rate. There¬ 
fore, sub-millimetre emission lines such as [CII] are an excellent 
way to observe the host galaxy in a frequency regime that may be 
less dominated by the SMBH. 

Using [CII] to investigate high-redshifQsources is a recently 
established field of study. iMaiolino et alj i2005l) were the first to 
obtain a [CII] detection in a high-redshift (z = 6.42) QSO, which 
they achieved using the IRAM 30-m single-dish telescope. Many 
early resolved detections of high-redshift [CII] were obtained 


with the IRAM Plateau de Bure Interferometer (PdBI; Walter et al. 

2009; 

Maiolino et al. 20121; iGallerani et al. 20121; IVenemans et al. 

2012; 

Cicone et al. 2015n. These observations have indicated that 


[CII] in high-redshift QSO host galaxies is found within 1 kpc of 
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the galaxy nucleus. The presence of [CII] on small spatial scales 
suggests that vigorous star formation is taking place in the cen¬ 
tral regions of these host galaxies, although a quantitative estimate 
of the overall star formation rate (SFR) based only on [CII] lu¬ 
minosity would be unreliable l lDIaz-Santos et aljl2013l) . IRAM ob¬ 
servations of [CII] revealed typical velocity widths of up to sev¬ 
eral hundred kilo metres per second. Ho weve r, in observat i ons of 
QSO J1148+5251. iMaiolino et all i2012h and lcicone et all fioTg) 
detected a second, much broader component of the line (FWHM ~ 
2000 km s -1 ) extending tens of kiloparsecs from the nucleus, sug¬ 
gesting the presence of a rapid, massive gas outflow powered by 
the QSO. 

ALMA is now poised to revolutionize [CII] studies of QSO 
host galaxies, as its incredible sensitivity allows high signal-to- 
noise detections with integration times of just a few minutes. In 
ALMA Early Science mode. lwillott et alld2013h detected [CII] at 
z = 6.4 in QSO CFHQS J0210—0456, obtaining similar results to 
the earlier IRAM observations. Willott et al. also obtained a new 
result: they observed a smooth velocity gradient of 100 km s -1 
across a spatial scale of 6 kpc, which c ould po te ntially originate 
from rotation of a galaxy-wide disc. Iwang et al] ( 120131) followed 
with a larger study in ALMA Cycle 0, detecting the [CII] line in 
five other z ~ 6 QSO host galaxies. They also observed a velocity 
gradient in four of their five targets, with line luminosities in the 
range 2-8 x 10 9 L 0 and FWHM in the range 200-500 km s _1 . 

These emerging high-redshift QSO studies will be crucial for 
probing the role that QSOs play in galaxy evolution. A correlation 
between the mass of the SMBH and the bulge stellar velocity dis- 
persion, known as the M-a relation (e.g.. lKormendv & Richstona 
119951: iGultekin et all 120091) . suggests a connection between the 
growth of a galaxy’s bulge and its SMBH, despite the fact that 
the radius of in fluence of a SMBH is sig nificantly smaller than 
the bulge radius dFerrarese & Merrittl feoOO). It is a fairly tight re¬ 
lation, b ut the amount and origin of the scatter remain controver- 
sial (e.g jGraham & Lil2009l : lGraham et alj201ll : lGraham & Scottl 
120151) . The most commonly proposed theory for this correlation is 
the presence of a ‘feedback’ mechanism, wherein the SMBH ex¬ 
erts radiative or mechanical energy, in proportion to its accretion 
rate, that slows or halts star formation in the host galaxy. The two 
most widely supported theories for the feedback mechanism are (a) 
radiatively driven gas outflows in luminous sources and (b) jet mo¬ 
mentum in radio-loud sources. An alternative theory to QSO feed¬ 
back is that the M-a relation arises naturally from the statistics of 
mass averaging in gal axy mergers 1 Pend 2007|; Ijahnke & Macciol 
l201lh . As presented in lKormendv & Hol d2013h . the local M-a re¬ 
lation is established from observations of about 50 nearby galaxies 
whose SMBH mass can be reliably determined through direct dy¬ 
namical measurements. Extension to high redshifts is tricky, as the 
relevant parameters must be determined through secondary meth¬ 
ods, and often involves an extrapolation that was originally cali¬ 
brated on the local relation. Where such high-redshift studies have 
been performed, they indicate that the ratio of SMBH mass to 
bulge properties is higher than the loc al ratio (e.g. JPeng et alj2006l : 
Wang et al. boidJ willott et al. l2013l) . Such results have been inter¬ 
preted as indication that SMBHs have an early, rapid growth phase 
that precedes galaxy grow th (e.g jMcLure et alj2006l : lReines et all 
1201 ll: I Willott et al.l2013l) . 

Studying the most luminous QSOs is another approach to 
these topics. Naturally, selection biases necessitate that high- 
redshift samples are intrinsically quite luminous; however, ex¬ 
tremely luminous QSOs are observed across a range of redshifts, 
although they are too rare to manifest at z < 1. The most lu¬ 


minous QSOs likely represent the most massive SMBHs under¬ 
going strong growth, and/or less luminous SMBHs experiencing 
super-Eddington accretion rates. While there is significant contro- 
versy as to whether QSOs are truly associated with mergers (e.g, , 
ITreister et al]|20ld : iKocevski et al.ll2012l : iHeckman & Bestll2014L 
and references therein), there is strong evidence that mergers do 
trigger the most luminous ones aTreister et al J 120121) . Addition¬ 
ally, there is evidence that strong gas outflows are ubiqui tous, or 
nearly so, in the host galaxies of lumin ous QSOs fe.g jGreen 2006; 
iKrongold et alj|20ld ; iLiu et alJl2014l) . Such outflows complicate 
these studies, as lines that indicate outflows are often the same 
emission lines that are used to estimate SMBH mass (e. g., CI V) 
or galaxy dynamics (e.g., [CIID at high reds hifts ( IRichards et al.l 
1201 ltlMaiolino et alj2012tlCicone et al.l2015t) . The most luminous 
sources are rare objects, pushing the boundaries of our theories and 
understanding, and thus are very interesting to study but also diffi¬ 
cult to interpret. 

In this manuscript, we present an ALMA Cycle 1 detec¬ 
tion of an extremely broad [CII] line in an extremely luminous 
(Lboi > 10 14 L 0 ) QSO. (Results and implications from the con¬ 
tinuum observations of our full QSO sample will be presented in a 
later publication.) Section [2] of this manuscript presents the target 
selection and observations. Section [3] presents the data reduction 
and analysis. In Section [4] we present interpretations of the line, 
provide estimates for the bulge and disc properties, and discuss this 
source in the context of the M-a relation. We conclude in Sec¬ 
tion [5] 

Throughout the remainder of this paper, we use a flat A cos¬ 
mology where Qa = 0.7 and Ho = 70 km s _1 Mpc -1 . 


2 TARGET SELECTION AND ALMA OBSERVATIONS 


We first identified luminous QSO s as p otential targets^ from the 
multi-wavelength QSO catalog of iKrawczvk et al] 1 20131) . which 
was compiled mainly from the Sloan Di gital Sky Survey (SDSS) 
Seventh Data Release Quasar Catalog I lSchneider et alJl2O10h and 
various broad-band sky surveys ranging from MIR to X-ray. Fol¬ 
lowing the recommendation of Krawczyk et al., we define a QSO’s 
bolometric luminosity as the energy in the (rest-frame) wavelength 
range 1 ^tm-0.62 nm (2 keV). This range, which is dominated 
by the accretion disc luminosity, corresponds to the intrinsic en¬ 
ergy of the QSO prior to re-processing of any kind, and there¬ 
fore gives the best estimate of bolometric luminosity when there 
is no significant covering fraction or absorbing material along the 
line of sight (as for Type 1 QSOs). Of ~ 120,000 optically identi¬ 
fied QSOs in the Krawczyk et al. catalog, fewer than 0.1 per cent 
have bolometric luminosity greater than 1O 14 L 0 , meaning that they 
are as or more luminous than the most luminous known obscured 
QSOs l lEisenhardt et alj|2012l; IWu et alj|2012l: iBridge et all 1 20131 
C. Lonsdale et al., in preparation). With a redshift of z = 4.612, the 
QSO SDSS J155426.16+193703.0 (hereafter J1554+1937) has the 
highest redshift of these most luminous QSOs from the Krawczyk 
et al. catalog. 

J1554+1937 has a bolometric luminosity Lbol = 1.3 x 
io 14 l q . There is no known evidence to suggest that J1554+1937 
is lensed (and thus has an artificially brightened apparent lumi¬ 
nosity). Other than its extremely high luminosity it does not seem 
to be an unusual QSO. Its broad-band spectral energy distribution 
(SED) and optical spectrum are typical for a high-redshift QSO. 
It is undetected in the NRAO-VLA Sky Survey dCondon et al.l 
11998h and in the Faint Images of the Radio Sky at 20 cm (FIRST) 
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Figure 1. We show 1900 GHz (158 fim) continuum (Panel A) and moment maps (as labelled; Panels B through D) of the [CII] line of J1554+1937. The grey 
dashed contours show the continuum profile in each panel (the color-scale and contours in Panel A are equivalent). The filled grey ellipse indicates the shape 
and orientation of the synthesized beam (analogous to ‘point-spread function’ in optical/IR observations). The emission shown in Panels C and D is restricted 
to regions where the signal-to-noise ratio in the [CII] data cube is greater than 4. 


survey (Becker et al.i 11995t) : thus its 1.4-GHz flux density is less 
than ~ 1 mjy. The typical criterion to classify a QSO as radio- 
loud or radio-quiet is the radio-to-optical ratio, R, originally de¬ 
fined in the observed frame: radio-loud QSOs generally have R 
val ues in the range 10-1000, while radio-quiet QSOs have R < 
1 (KellermannetaPl 19891) . J1554+1937 has an observed-frame 
radio-to-optical ratio R < 40 (determined from the FIRST sur¬ 
vey upper limit and the g -band magnitude), meaning that it may be 
radio-quiet, radio-intermediate, or at the low end of the radio-loud 
range. 

The redshift of J1554+1937 makes its [CII] line accessible 
to ALMA’s Band 7 (275-373 GHz). QSO emission lines are quite 
broad (ty pi cally > 1000 km s _1 even for the narrow lines; e.g., 
iHao et al.ll2005r) . and therefore do not require high spectral resolu¬ 
tion. We placed the centres of two of the 2-GHz spectral windows 
at 337 and 339 GHz, covering the full 4-GHz-wide lower sideband 
with a spectral resolution of ~ 30 km s _1 . A small frequency range 
(19.75 MHz) at 338 GHz is unusable due to low sensitivity at the 
spectral window edges. Two other spectral windows covered the 
frequency range 348-352 GHz. 

In interferometric observations, there is a necessary trade-off 
between better angular resolution (with a more extended array con¬ 
figuration) and better sensitivity (with a more compact array config¬ 
uration). The primary motivation for our observations was a contin¬ 
uum detection experiment (the results of which will be discussed in 
a forthcoming paper), and for this reason we opted for an array con¬ 
figuration that would achieve high sensitivity at the cost of spatial 


resolution. The observations were performed using a compact con¬ 
figuration, resulting in angular resolution of approximately 1 arcsec 
at 340 GHz. 

J1554+1937 was observed for 5.2 minutes with 25 available 
antennas. Titan was used as the absolute flux calibration source and 
nearby QSO J1619+2247 (6.6 degrees distant) as the phase calibra¬ 
tor. 


3 DATA ANALYSIS 

Initial data reduction was performed by staff at the North American 
ALMA Science Cente{]] (NAASC). The data package provided to 
us by the NAASC included Python scripts that would execute the 
basic data reduction and calibration; the data package also included 
preliminary images^The scripts were executed in version 4.2 of 
CASA0 (McMullin et al] 120071 ). We performed minor modifica¬ 
tions to the scripts in order to reduce the number of frequency chan¬ 
nels flagged at the spectral window edges. We produced our own 
set of de-convolved images using CASA’s clean task, for which 
we used Briggs weighting with parameter setting robust = 0.5. 
We performed a single round of phase-only self-calibration using 
solutions from the combined continuum and line data. The synthe¬ 
sized beam size (analogous to the ‘point-spread function" in opti- 

1 https://science.nrao.edu/facilities/alma 

2 http://casa.nrao.edu/ 
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cal/IR observations) is I”2 x 0"7, corresponding to 7.8 x 4.6 kpc 
at 2 = 4.6. 

The continuum map (in the frequency ranges 336-337 GHz, 
339-340 GHz, and 348-352 GHz) is shown in Panel A of Fig.Q 
the rms noise level in the continuum image is approximately 
0.15 mJy beam -1 . To estimate the source size, we performed a 
Gaussian fit to the data in the UV-plane using the CASA task 
uvmodelf it, which yielded a de-convolved source size of 310± 
30 x 260 ± 70 mas (~ 2.0 x 1.7 kpc) indicating that the source 
is marginally resolved spatially. Based on the maximum baseline 
length of ~ 400 m, the smallest angular scale we could resolve at 
338 GHz is O'/46. The integrated flux density resulting from the fit 
is 12.70 ± 0.12 mJy at ~ 344 GHz. At a redshift of z = 4.6, that 
flux density corresponds to a luminosity of ~ 2.3 x 10 12 Lq at a 
rest frequency of 1930 GHz (156 /im). This value is three orders 
of magnitude more luminous than the upper limit of the 1.4-GHz 
luminosity (see Sec. [2}; therefore it is clear that the 1930-GHz con¬ 
tinuum emission is due to dust rather than a radio jet. 

We isolated the line emission by subtracting the continuum 
from the data cube. An intensity map (moment 0 map) of the line 
using the frequency range 337.3-338.2 GHz is shown in Panel B 
of Fig. [I] the rms noise level in this image is approximately 
0.4 Jy beam -1 km s -1 . We again used CASA's uvmodelfit 
task to estimate the spatial extent of the line-emitting region. 
Formally, the line image is resolved with de-convolved size of 
250 ± 40 x 250 ± 140 mas (1.6 x 1.6 kpc), which is consistent 
with the size estimate of the continuum emission. The integrated 
flux density is to be 9.40 ± 0.40 Jy km s -1 , corresponding to a line 
luminosity of 5.9 x 10 9 Lq. 

The velocity vs. position map (moment 1 map) of the [CII] 
line is shown in Panel C of Fig. [I] We have restricted this analysis 
to pixels in the data cube where the S/N is at least 4x the rms noise 
level of ~ 2 mJy (determined in a line-free channel). The velocity 
map shows a clear gradient of higher to lower velocity from the 
east to the west. In the moment 0 map, the S/N is ~ 10 at the edges 
of this region, indicating that this gradient is real. This velocity 
gradient is suggestive of large-scale, ordered rotation; implications 
are discussed in Section[4] 

In Panel D of Fig. |T| we present the velocity dispersion map 
(moment 2 map) of the [CII] line. The central emission region has 
a velocity dispersion of about 220 km s -1 . On the edges of the 
emission region, the velocity dispersion drops to < 100 km s -1 . 
In Section |4~71 we discuss the meaning of this velocity dispersion 
according to whether it originates from the host galaxy bulge or 
from a rotating star-formation disc. 

The [CII] line profile is shown in Fig. [2] The expected red- 
shift was 2 = 4.612 (based on the optical spectrum). The actual 
redshift of the line is approximately 4.627: the observed line centre 
is 337.8 GHz (887.6 gm). Because of this offset from the expected 
frequency, the line falls over the edges of two of ALMA’s spectral 
windows; approximately 19.75 MHz (6.7 km s -1 ) of data are miss¬ 
ing (flagged) at 338 GHz. The profile is clearly non-Gaussian. It is 
closer to being flat-topped but also shows line structure; it may be 
a double-horned profile. The FWHM is 735 km s -1 . 

4 ANALYSIS 

4.1 Redshift offset of the C IV and [C II] lines 

We find a redshift offset between the [CII] emission line and the 
UV lines in the SDSS spectrum. The [CII] line is at 2 = 4.627 


observed frequency [GHz] 


340 339 338 337 336 



relative velocity [km/s] 


Figure 2. ALMA detection of [CII] line in J1554+1937. The rms varia¬ 
tion (determined from line-free channels) is approximately 2 mJy. Observed 
frequencies are indicated on the top axis. The lower axis shows the velocity 
relative to the line centre based on a redshift of z = 4.627. The dashed 
vertical line indicates the line centre; the dotted vertical line indicates the 
expected line centre for a redshift of z = 4.612 (the visually confirmed 
SDSS redshift). Missing data at 338 GHz correspond to the low sensitivity 
edges of the spectral windows. 


(based on the central frequency of 337.8 GHz), while CIV, NV, 
and Ly a a ppear to be at lower redshift. The SDSS pipeline 
dBolton et alj2012h measured the redshift of Ly a to be 2 = 4.622; 
however, its true redshift is difficult to pinpoint owing to asym¬ 
metric absorption by the Ly a forest. We consider the ‘true’ red¬ 
shift of the ultraviolet lines to be the visually determined redshift 
reported in SDSS Data Release 9 lAhn et alil2012l ). The redshift 
determined by eye is 2 = 4.612, corresponding to a velocity off¬ 
set of ~ 800 k m s -1 from [CIIl, Su ch o ffsets are not unusua l 
in QSOs. Indeed. [Sulentic et afl (2000 ) and iRichards et all J2002h 
have observed that the CIV line is nearly always blue-shifted with 
respect to low-ionisation lines, with offsets that occasionally ex¬ 
ceed 2000 km s -1 , although 800 km s -1 is more typical. 

One theory put forth to explain the large blueshift in broad, 
high-ionisati on lines like CIV is commonly referred to as the ‘disc- 
wind’ model dSulentic et al.l200ol ; lRichards et alj201 lh . According 
to this picture, the line is emitted from a radial outflow whose far 
side is obscured by the accretion disc, and thus only the blue-shifted 
component is visible. The wind may arise from the disc itself. 

We wish to emphasize this caveat to future ALMA observers: 
in the sub-millimetre regime, such velocity offsets can amount to 
shifts of ~ 1 GHz in the observed frame, potentially moving a tar¬ 
geted sub-millimetre line halfway across one of ALMA’s 2-GHz- 
wide spectral windows. If the redshift of a QSO is known only from 
the (rest-frame) UV lines, optical emission lines, we therefore rec¬ 
ommend that ALMA observers place multiple spectral windows 
alongside each other to allow for an unknown, but potentially quite 
large, shift in the line frequency. 


4.2 Implications from the [CII] velocity map 

The velocity map of the [CII] line (Panel C of Fig. [TJ shows a gra¬ 
dient across the source, with the eastern side of the source being 
redshifted with respect to the western side. This morphology sug¬ 
gests large-scale rotation of a gravitationally bound gas component. 
Similar line morp holo gies were observed i n five 2 > 6 QSOs by 
IWang et all d20 1 3l) and Iwillott et al.l d2013h . In those cases, as for 
the source presented here, the observations were not high-spatial- 
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resolution, such that the detailed velocity morphology is not tightly 
constrained. However, the overall gradient is clear. 

One possible origin for such a velocity gradient is a rotat¬ 
ing gas disc that is observed partially edge-on. Is the velocity dis¬ 
persion map (Panel D of Fig. Q} consistent with a rotating disc? 
Generally, a rotation-supported disc should show a velocity disper¬ 
sion that is fairly constant across the disc, representing the intrin¬ 
sic dispersion of the disc material. For an inclined disc, however, 
poorly resolved rotation can cont aminate the dispersion map, an ef¬ 
fect k nown as ‘beam-smearing’ toavies et all201 ll : I New man et al.1 
120131) . A velocity-dispersion peak near the centroid of an unre¬ 
solved source strengthens the evidence for rotating gas. This mor¬ 
phology is in fact what we observe for J1554+1937 (Panel D of 
Fig.Q}. which has a velocity dispersion peak of ~ 220 km s 1 near 
the source centre. The true intrinsic dispersion is better measured 
near the outskirts of the emission region, which in this case suggest 
an intrinsic dispersion of 50-100 km s -1 if the emission is truly 
from a disc. This level of dispersion is similar to val ues measured 
for mo lecular gas discs by Swin bank et alj J20I ll) and lHodee et al.l 
l l2012l) . However, there are problems with such an interpretation. 
Studies of high-redshift galaxy morphology with the Hubble Space 
Telescope argue against the p resence of well ordered galaxy discs 
at redshifts above ~ 2 (e. g., iGlazebrook et al.ll 19951 : iDriver et al.l 
ll998l ; lFerguson et al.l2004l :TPaoovich et al.l2005l) . At high redshift, 
star formation seems to have occurred in massive clumps (lC- 
10 9 Mo ) within turbulent discs dElmegreen & Elmeg ree nl [20051: 
Elmeg reen et al ] [2009 ). Theoretical formulations (e.g., Dekel et all 


2009) support the notion that stable discs did not form in galaxies 
before 2 ~ 1 . 

If the [CII] emission does not originate in a rotating disc, it 
may instead be coming from the host galaxy’s bulge, or a prog enitor 
to the bulges that are seen in local galaxies (e.g.. lElmegreen et al.l 
120091) . which would also explain the small spatial extent of high- 
redshift [CII] emission. Galaxy bulges in the local Universe are 
typica lly ‘red and dead’ (i.e., with little or no gas or star-forma tion) 
(e.g.. lRenzinill2006l : lBlanton & Moustakasll2009l : lFabianl2012h . but 
at higher redshift they must have undergone intense star formation 
in order to generate the massive stellar bulges observed today. The 
luminous [CII] line in this and other sources may be due to vig¬ 
orous star formation in the young bulges, or bulge progenitors, of 
early-stage galaxies. 

Bi-directional outflows of ionized gas are another possible ori¬ 
gin for the [CII] velocity gradient. Outflows have been claimed 
to be ubiquitous in luminous QSO s, as evidenced by mapping of 
extended ionized-gas regions (e.g.. lLiu~et ai1 l20l4 McElrov et al. 


2015) and anal ysis of broad absor ption lines (e.g., IShen et al 


2008 ). In fact, iGangulv et al.l d2007l) observed a positive correla¬ 


tion between the presence of broad absorption lines and a QSO’s 
luminosity, suggesting that outflows are m ore pr evalent in the 
most luminous QS Os. Furthermore, iMaiolino et al] J2012h and 
ICicone et alJ <1 20151) observed a striking [CII] line profile in a 
2 = 6.4 QSO (one of the most distant known, with bolometric 
luminosity 8 x 10 12 Lq), with broad wings that are attributed 
to a massive outflow. While the narrow component of the [CII] 
line in that source has FWHM of ~ 350 km s -1 , the low-level 
broad wings are consistent with a Gaussian profile extending to 
±1300 km s“\ and are marginally resolved on the physical scale 
of the host galaxy (~ 16 kpc). If the [CII] emission in J1554+1937 
originates from outflowing material, that material must have a large 
range of velocities in order to produce the observed line profile 
shown in Fig. [2] However, as the profile and velocity gradient can 
also be explained by the presence of a massive rotating disc, in- 
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yoking outflows as a possible explanation is not strictly necessary. 
Perhaps the [CII] emission is actually generated by a combination 
of disc rotation and outflowing gas. 

Another physical process that is associated with luminous 
QSOs is that of major mergers between gas-rich galaxies. While 
QSO activity can be induced by secular (internal) processes within 
a galaxy, major mergers seem to be the driving force behind 
the most luminous QSOs l lSanders & Mirabel 1 19961 : fUrrutia et al.l 
1 20081 : FlVeister et all2012l) . If a merger was the precursor to this ex¬ 
tremely luminous QSO, then the gas kinematics likely show signs 
of disruption or turbulence, which could be revealed by higher- 
spatial-resolution spectroscopy. Activity resulting from a merger 
does not remove the possibility that outflowing gas is present. On 
the contrary, the evolutionary framework describing QSO forma¬ 
tion purports that major galaxy mergers trigger a ‘feedback’ phase 
where energy fr om the SMBH disrupts the host galaxy’s star for- 
mation activity (Sanders et al. 1988; Kauffmann & Haehnelt 2000; 


IWvithe & Loebll2003l : Di Matteo et alj|2005l ; Hopkins e t al. 2008); 

outflows are one possible form of QSO feedback. 

The obvious next step toward identifying the true nature of 
the [C II] line is to obtain higher-spatial-resolution spectroscopy 
with ALMA. Such observations would allow the isolation of the 
line profile in different ph ysica l regions of the galaxy, a technique 
known as ‘kinemetry’ (e.g. lKrainovic et ap|2006i) . Asymmetries in 
such mappings have been used to di stinguish between disc rotation 
and merger activity ( Shapiro et al. 2008), while ioni sation diagnos- 
tics from kinemetry can trace directional outflows IMcElrov et al J 
[20151) . 


4.3 [CII] line profile 

Based on the velocity map, we have speculated that the line orig¬ 
inates in a massive rotating disc. In this section, we show that a 
rotating disc can generate a line profile similar to the one shown in 
Fig. 0 We avoid complex modelling at present, as the low spatial 
resolution of the current observations does not constrain the possi¬ 
ble dynamical theories. 

In general, the line profile from a rotating disc depends on 
the velocity field (i.e., the rotation curve), the surface brightness 
profile (i.e., the gas distribution), and the inclination of the disc. 
The surface brightness profile and the inclination determine the line 
intensity at each position on the projected disc image, while the 
rotation curve determines the observed velocity of the emission at 
each position. 

Observed emission lines from discs show a h uge range of pro- 
files, as demonstrated by neutral hydrogen (HI) ( IKoribalski et al.l 
2004). The two extreme types of line profile are the singly peaked 
(typically Gaussian) profile of irregular dwarf galaxies and the 
double-peaked or ‘double-horned’ profiles of large spiral galax¬ 
ies. Both types of galaxy have exponential surface brightness pro¬ 
files, but their velocity fields are different. The rotation curve of a 
typical spiral galaxy initially i ncreases with radius, then levels of f 
at some limiting velocity (e.g.. lRubin et alJ[l978l : lcremen 3 11983) . 
Such galaxies have a double-horned HI profile due to a large range 
of radii generating emission in a small range of velocities, result¬ 
ing in symmetric red-shifted and blue-shifted peaks corresponding 
to the maximum rotational velocity offset from the systemic ve¬ 
locity. In contrast, late-type dwarf galaxies typically have a con¬ 
tinually rising velocity curve (e.g. JCote et alj|2000t) . resulting in a 
single-peaked line at the systemic velocity. In between those two 
extremes, it is possible to have a nearly flat-topped profile. This 
profile manifests when the characteristic radius (the ‘knee’) of the 
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v - v sys (km/s) 

Figure 3. Examples of three different line profiles from dynamically equiv¬ 
alent discs. The formalisation of the model is given in the text. For all three 
profiles here, t>o = 480 km s _1 , rr v = 75 km s —1 , and i = 60°. The 
value r c /ro = 0.5 generates the single-peaked (red) profile; r c /ro = 1 
generates the flat-topped (brown) profile; r c /ro = 2 generates the double- 
peaked (black) profile. The line amplitudes are arbitrary for the purposes of 
this figure. 


surface brightness profile is the same as the characteristic radius of 
the velocity rotation curve. In other words, the surface brightness 
falls roughly as the velocity curve rises, resulting in constant line 
luminosity across a range of observed velocities. 

We demonstrate idealized line profiles using a disc-galaxy 
model with standard rotation curve and surface brightness distri¬ 
bution, as follows. We assume a fiducial velocity (v) rotation curve 
that initially rises exponentially, then levels off with some charac¬ 
teristic radius ro such that v(r) = «o[l — exp(—r/ro)]. For the 
surface brightness (/), we use an exponential profile with charac¬ 
teristic radius r c such that 7(r) oc exp(— r/r c ), as is seen in nor¬ 
mal disc galaxies. For perfectly ordered rotation, the emission line 
from any point on the disc could be approximated as a delta func¬ 
tion. However, the addition of an intrinsic velocity dispersion (a v ) 
yields a line profile ( L ) at any point on the disc given by 


L(v) 



e 




2 


( 1 ) 


where v(r ) corresponds to the average velocity at that radius. Fi¬ 
nally, we allow for an inclination of the disc to the line of sight. 
The spatially integrated line profile from the disc is therefore the 
convolution of the line profile and the brightness profile integrated 
over the (projected) face of the disc. 

Example line profiles are shown in Fig. [3] For the purposes of 
this demonstration, we assign the parameters vo = 480 km s -1 , 
a v = 75 km s _1 , and i = 60°. A ratio of r c /ro = 0.5 gener¬ 
ates the single-peaked profile, r c /ro = 1 generates the flat-topped 
profile, and r c /ro = 2 generates the double-peaked profile. These 
three profiles thus correspond to discs with equivalent dynamics; 
the variation is due to the ratio of r c /ro. In other words, the ob¬ 
served profile depends on the surface brightness at the outer edge 
of the disc, where the rotational velocity is at its maximum. 

A ratio of r c /ro > 1 is required but not sufficient to gener¬ 
ate a double-peaked line profile. A single-peaked profile can result 
from an intrinsic velocity dispersion that is high enough to mask 
the double-horned rotation signature. This effect is demonstrated in 
Kg. a which shows three different models corresponding to vary¬ 
ing values of maximum velocity (vo) and intrinsic velocity disper¬ 
sion (<r„). Relatively higher intrinsic velocity can dominate over the 
wings of the line associated with rotation, yielding a Gaussian pro¬ 



Figure 4. Line profiles shown here demonstrate that a large intrinsic ve¬ 
locity dispersion in the disc can mask the double-horned rotation signa¬ 
ture, resulting in a single-peaked profile. The formalisation of the model 
is given in the text. For the three profiles here, r c /ro = 2 and i = 60°. 
The maximum velocity (uo) and intrinsic velocity dispersion (a v ) in kilo¬ 
metres per second are: no = 350, a v = 175 for the single-peaked (red) 
profile; vo = 415, cr v = 110 for the nearly fiat-topped (brown) profile; 
and vo = 480, <r v = 75 for the strongly double-peaked (black) profile. 
The purple data points correspond to the [CII] line profile of J1554+1937 
shown in Fig. [2] 


file. The fact that the [CII] line observed in J1554+1937 appears 
to be double-peaked or flat-topped rules out the notion that the line 
width is merely a result of high velocity dispersion. 

It is important to emphasize that the parameters for inclination 
(i) and maximum velocity (vo) are degenerate in this model. There¬ 
fore, although the model formulated above separates these two pa¬ 
rameters, we do not claim to know the true parameters as a result of 
these rough fits. Qualitatively, the profile shape implies that, if the 
line is generated by a disc, then the disc has a relatively high maxi¬ 
mum velocity and a relatively low intrinsic velocity dispersion. 


4.4 Estimating dynamical mass for a rotating disc 

If we assume a rotating disc geometry for the origin of the [CII] 
line, we can estimate the dynamical mass within the disc radius 
as Moyn / Mq ~ 1.16 x 10 5 Ho D where vo is the maximum cir¬ 
cular velocity of the disc in km s -1 and D is the disc diam¬ 
eter in kpc. We visually estimate a maximum circular velocity 
using the half-width of the line at 20 per cent of its maximum 
value: vo x sin(i) ~ 480 ± 50 km s , where i is the incli¬ 
nation angle of the disc relative to the line of sight (i = 0° for 
a face-on disc). Using the source size from the line data fitting 
{D ~ 1.6 kpc; see Section [3). the estimated dynamical mass is 
Mdyn x sin 2 (i) ~ (4.3± 1.2) x 10 10 Mq. Using the size estimate 
from the continuum data (D ~ 2.0 kpc) yields a slightly larger 
mass estimate of Ma yn x sin 2 (i) ~ (5.3 ± 1.3) x 10 10 Mq. The 
poor spatial resolution prevents us from determining the inclina¬ 
tion, and shape-based inclination estimates are controversial in any 
case due to unknown intrinsic ellipticity and variations of elliptic- 
ity with radius jBellovarv et alJl2014l) . Certainly, unless the disc is 
very close to fully face-on (e.g., i < 20°), the dynamical mass must 
be in the range of 5 x 10 10 -6 x 10 11 Mq. Thus, this source prob¬ 
ably has a d ynami cal mass similar to the a ~ 6 QSOs studied by 
IWang et al] (12013 1); suc h values are also typical of QSOs at lower 
redshifts ( iSolomon & Vanden Boulli2005l) . 

We can compare the estimated dynamical mass to that same 
parameter observed in other gas-rich galaxies, both locally and 
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at high redshift. Rotation curves and dynamical masses are of¬ 
ten determined using hydrogen lines such as Ha, which probes 
the star-forming regions of a galaxy. For example, the Sydney- 
AAO Multi-object Integral field spectrograph (SAMI) Galaxy Sur¬ 
vey ( lAllen et al J|20 15T) has produced high-angular-resolution data 
cubes for a large number of local galaxies with a range of galaxy 
types. In about 120 galaxies from this sample with z < 0.06, 
G. Cecil et al. (in prep.) finds maximum dynamical masses of about 
10 11 M w , with a typical value of about 2 x 10 10 M 0 . At 2 ~ 2, 
IForster Schreiber et al .1 1200a) observed Ha in several UV-selected 
galaxies and found that they were generally consistent with having 
clumpy, rotating gas discs within the half-light radius of ~ 4 kpc, 
with dynamical masses ranging from ~ 0.5 to 25 x 10 10 Mq . Thus, 
the mass we estimate for 11554+1937 places it in the high-mass tail 
of the galaxy population. 

4.5 Estimating central velocity dispersion for a bulge 

If the [CII] line originates from the host galaxy bulge, then we can 
use the velocity dispersion of the line to investigate the source with 
respect to the M-a relation: the empirical scaling relation that re¬ 
lates the mass (A/bh) of a galaxy’s SMBH to the stellar velocity 
dispersion (<r) of its bulge component. We first estimate the bulge 
velocity dispersion in this section; the results are discussed in con¬ 
text of the M-a relation in Section l4~7l 

To estimate the stellar velocity dispersion, the simplest ap¬ 
proach is to use the velocity dispersion of the [C_II] li ne di¬ 
rect ly, a method that w as also used by IVenemans et al J 1 201 2ll 
and lWillott et al J (1 201 3l) for high-redshift QSOs detected in [CII], 
For J1554+1937, the velocity dispersion map suggests that a ~ 
220 km s -1 (see Panel D of Fig. G). We caution that the value 
we measure may not be the intrinsic bulge velocity dispersion if 
there are any non-bulge co mponents contributing t o the measured 
line dispersion. For example. lBellovarv et al.l 02O144 show how disc 
rotation and inclination effects combine to falsely inflate disper¬ 
sion measurements. Additionally, any rotating component that is 
not spatially resolved can incre ase the velocity dispersion measure¬ 
ment due to ‘beam-smearing’ jDavies et alJuOl ll : iNewman et al.l 
120131) . These effects are extremely difficult to account for without 
knowing the true galaxy inclination. 

Does the dispersion of the [CII] line provide a reliable mea¬ 
surement of the bulge stellar velocity dispersion? For most of the 
other distant sources with estimated stellar velocity dispersion, that 
value was deter mined from the widths of narrow UV/optical emis¬ 
sion lines (e.g.. IShields et all 120031 : ISalviander et alj [20071) . The 
[OIII] forbidden line, for example, originates from the narrow¬ 
line region surrounding the black hole, and its width correlates 
well with the stella r velocity dispersion in the bulge dNelsonl 
l200d B orosonl 12003). That correlation suggests that the velocity 
field of the narrow-line region is determined in large part by the 
same gravitational dyn amics that govern the stellar velocity field 
dNelson & Whittle!! 1 9961) . The question of whether [CII] similarly 
probes the stellar velocity dispersion depends on the nature of the 
[C II]-emitting gas. As discussed above, it is likely that a young 
galaxy bulge at z = 4.6 will be experiencing significant star- 
formation. Th e [CII] luminosity is correlated with SFR in norma l 
galaxies (e.g.. iHollenbach & Tielensl 1 1 9971 : Ide Looze et al,ll201lh . 
but in active galaxies its luminosity is thought to be ‘contaminated" 
by g as that has been ionized by radiation from th e active nucleus 
(e.g jMalhotra et alJl997l ; lDfaz-Santos et alj2013l) . The [CII] line 
need not sample the bulge at all if the bulge region is highly ionized. 

We now return to the possibility that [CII] originates from a 


rotating star-forming disc instead of from a bulge region, as dis¬ 
cussed in Sections 14.21 and 14.31 If the [CII] line traces a rotating 
disc, then the velocity dispersion map is likely a result of ‘beam- 
smearing’. Even so, it may be possible to obtain a velocity disper¬ 
sion estimate f or the bulge component u sing the following relation. 
As discussed in lNelson & Whittle! dl 9961 ). the velocity dispersion of 
the bulge measures its gravitational potential while the maximum 
disc rotation velocity measures the potential of the disc and/or halo. 
Thus, simple theory based on the mass density profiles of galaxies 
predicts that the ratio between the maximum circular velocity of a 
galaxy’s disc and the velocity dispersion of its bulge should fall in 
the range 2 1 / 2 -3 1 / 2 , with the exact value depending on the mass 
profile. 

IFerraresei d2002l> was the first to empirically quantify the rela¬ 
tionship between circular rotation velocity (vo) and central veloc¬ 
ity dispersion (a). Using a nearby sample of 33 elliptical and spi¬ 
ral galaxies, Ferrarese derived a power-law fit in the region where 
70 < a < 180 km s -1 , but noted that for elliptical galaxies it may 
also be valid in the r ange 1 20 < a < 350 km s -1 . Following that 
result. lPizzella et al] 120051) examined a broader sample of 72 ellip¬ 
tical and spiral galaxies to investigate the relation between a and 
Vo separately for high- and low-surface-brightness galaxies. They 
found that these two classes of galaxies follow different relations, 
although the latter was determined from a sample of just eight. 
In addition, Pizzella et al. found that the high-surface-brightness 
and elliptical galaxies follow a single relation that is better charac¬ 
terised by a linear fit in the vo -a plane than by the power-law fit of 
IFerraresei 02OO2I) . Their linear relation fits their data well for veloc¬ 
ity dispersion as low as a ~ 50 km s^ 1 , and is a much better fit to 
the data with 300 < a < 350 km s -1 (cf. their Figure 2) than the 
relation fro mlFerraresel d2002l) . We therefore use the linear relation 
from lPizzella et alj d2005h ~to represent the relationship between a 
and vo- We reproduce their linear fit here as 

( 7 -^) = (1.32 ± 0.09) (-^1-)+(46 ±14), (2) 

V km s / V km s / 

where vo is determined where the rotation velocity curve is flat, 
and a e / s signifies that the value of a corresponds to an aperture of 
size r e / 8 . The radius r e is the “effective radius” of the bulge — the 
radius that contains half of the total amount of light dBaggett et alj 
Il998h . 

Without knowing the actual inclination of the disc, we can¬ 
not determine the true de-projected maximum rotation velocity, vo. 
However, we can estimate the maximum projected circular velocity 
from the width of the [CII] line. A reasonable estimate for the pro¬ 
jected maximum rotational velocity is the half-width of the line at 
20 per cent of its maximum, yielding vo = 480 km s -1 / sin(i); a 
profile corresponding to this value is shown in Fig.[4]as the double¬ 
horned black line. Applying this value in Equation[2]yields the ap¬ 
proximation that a e / s as 330 km s _ 1 /sin(i). However, the pro¬ 
jected maximum rotation velocity is a lower limit to the true max¬ 
imum rotation velocity, such that vo > 480 km s -1 ; similarly it 
follows that a e /g > 330 km s -1 . 

In Section 14771 we compare the prope rties of J 15 54+1937 to 
the local M-a relation as presented by iKormendv & Hoi J2013t) . 
Those authors use a value of a corresponding to an aperture of size 
r e / 2 , and therefore we must perform an aperture correction to_our 
a e / 8 estimate. We apply the correction given by Ijorgensen et all 
d 19951) . who used kinematic models based on empirical data to 
derive a formula for the measured central velocity dispersion as 
a function of aperture size. Using their formulation, we find that 
&e /2 = 0 . 93 x< 7 e / 8 ; therefore <j e /2 > 310 km s _1 forI1554+1937. 
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For inclination i = 60°, we find <r e /2 ~ 360 km s^ 1 , and for 
i = 45°, we find a e / 2 ~ 450 km s _1 . 

The largest known galaxy bu lge velocity dispersion s are < 
400 km s -1 (see references within IKormendv & Holl20130 . which 
means that J1554+1937 source is in the high-velocity-dispersion 
tail of the known distribution. A virialized bulge with a ~ 
400 km s -1 must have a high mass or a small radius, or both. For a 
simple estimate of the properties of such a bulge, we use the virial 
theorem, which states: 


Aft 


bulge ' 


3 Ra 2 
G 


(3) 


where Mb u i ge is the bulge mass and R is the virial radius. For 
this system, the bulge mass must be less than the total dynamical 
mass, which was calculated to be Md yn x sin 2 (i) w (5.3 ± 1.3) x 
10 10 Mg. If the bulge has virialized, then Equation [3] leads to the 
conclusion that R < 0.8 kpc. Such a bulge is within the typical 
population distribution for galaxies in the redshift range 1 < 2 < 3 
IBruce et al 1 1120141 ). although toward the high-mass/low-radius side 
of that distribution. Such compact cores are rare in nearby galax¬ 
ies, but it has been demonstrated that the cores of galaxies at high 
redshift are usually mor e compact than their local counterparts 
dvan DokkumetaT 20141). 

From IKormendv & Hoi d2013h , we note an im portant caveat 
to this analysis. In contrast to the iFerraresel { 20021 ) results, Kor- 
mendy & Ho assert that the tight correlation between vo and a 
only exists in the range of 150 km s^ 1 < Vo < 300 km s -1 . 
At higher rotational velocities, data are scarce and the scatter visi¬ 
bly increases. Therefore, our estimate for bulge velocity dispersion 
must be taken cautiously. Furthermore, the correlation holds for 
galaxies with classical bulges (having the light profile of an ellipti¬ 
cal galaxy, resulting from galaxy mergers) but not galaxies with so- 
called “pseudobulges” (resulting from secular processes). Applying 
the vo /<t relation to J1554+1937 thus requires the assumption that 
its host galaxy contains a classical bulge. 


4.6 Measuring the central black hole mass 

IShen et ail 1201 ll) estimated the SMBH in J1554+1937 to have a 
mass of Mbh = (2.8 ± 0.5) x 10 10 Mg. This mass estimate was 
based on the the w idth of the C IV line according to the calibration 
presented by IVestergaard & Petersonl ( 120061) . The error published 
by Shen et al. represents only the measurement error from the 
line fitting procedure, which equates to an uncertainty of 0.08 dex. 
However, Vestergaard & Peterson assert that the uncertainty in their 
mass calibrations based on CIV width is actually 0.32 dex. Based 
on the latter estimate for uncertainty, we consider the probable mass 
range for the SMBH in J1554+1937 to be (1.3-5.9)x 10 10 Mq. 

tPark et alj I 2 OI 3 I) presented an updated mass calibration for 
CIV that would decrease the SMBH mass e s timate to Mbh = 
6.3 x 10 9 Mq. However. iKratzer & Richards] ( i2015t) showed that 
this new calibration is likely to be systematically low by 0.5 dex 
or more, a conclusion that supports the original mass estimate for 
this source. While such a mass is toward the high end of what 
is theoretically observable dKellv et alJhoiol) . it is not unprece¬ 
dented, as other black holes approximately thi s mass ive are known 
at high redshift (based on other CIV widths; IShen et ai]|201ll) as 
we ll as locally (based on small-scale stellar velocity dispersions; 
iMcConnell etai]|201ll) . 

Does the CIV line width provide a reliable measurement of 
the black hole mass ? This line is thought to be the least reliable 
estimator (e.g.. iHo et al.|]20O ). but at high redshifts it is usually 


the only estimator available. For nearby (< 100 Mpc) systems, 
the SMBH mass can be reliably determined by dynamical mea¬ 
surements, whereas higher-redshift characterisations are based on 
single-epoch measurements of broad-line widths (e.g., Ha, Mgll, 
or CIV) that were first calibrated against the local M-a relation. 
Owing to a selection bias at high redshift, single-epoch esti mato rs 
may falsely conflate SMBH mass estimates ( IShen & KellviuOlOl ). 
Additionally, the CIV line width might be influenced by outflows, 
particularly in sources with large offsets between the CIV red¬ 
shift and the systemic redshift (see Sec. l4.lt . or by a major galaxy 
merger. On the other hand, iFine et al.l d20ld) have argued that 
single-epoch estimators actually depend more on intrinsic QSO lu¬ 
minosity than on the line widths themselves. In this case, despite 
the inherent uncertainty in using CIV as a mass estimator, the 
high bolometric luminosity of J1554+1937 indicates that its SMBH 
must be extremely massive. If the SMBH is Eddington-limited (as¬ 
suming the classical optically thin case), then its mass is actually 
much lower: ~ 4.0 x 10 9 Mq. However, such a luminous system 
could also result from an even lower-mass SMBH that is radiat¬ 
ing at a super-Eddington rate, thoug h su ch super-Eddington QSOs 
are expected to be quite rare dKellv et al |20 1 (II . If the current mass 
estimate is correct, then the Eddington ratio (-Lboi/T/Edd) is 0.15. 


4.7 Comparison to the M-a relation 

With estimates in hand for the bulge velocity dispersion (Sec. 1431 
and central black hole mass (Sec l46t , we are now equipped to 
compare J1554+1937 to the local M-a relation dKormendv & Hoi 
|201 3l) . which we reproduce here as 


log 


f Mbh \ 

\M^J ~ 

(8.49 ± 0.05) + (4.38 ± 0.29) log i) ' (4) 


In Figure[5] we show the data and fit presented by Kormendy & Ho. 
We show the estimated locations for J1554+1937 for different po¬ 
tential values of the unknown disc inclination, using the black hole 
mass estimate of Mbh = (2.8)( 9 'g) x 10 1 o Mq. For the minimum 
estimated value of a e /2 = 310 km s -1 (corresponding to i = 90°), 
we see that J1554+1937 would be over an order of magnitude more 
massive in Mbh than predicted by the local relation. Given that the 
intrinsic scatter of the M-a relation is 0.29 dex dKormendv & Hoi 
|201 3l) . such an offset corresponds to a discrepancy at the level of 
2.6a rms . For inclinations in the range 27° < i < 45°, J1554+1937 
would consistent with the local M-a relation within l<r rma . There¬ 
fore, for likely inclinations between 45° and 90°, J1554+1937 can 
be considered over-massive with respect to the relation by an order 
of magnitude or less, but not at a significant level. 

For comparison, we also show in Figure [5] the location of 
J1554+1937 if one were to equate the velocity dispersion of the 
[CII] line directly with the velocity dispersion of the bulge. In tak¬ 
ing that approach, we would find that J1554+1937 is nearly two 
orders of magnitude above the local M-a relation (with a signif¬ 
icance of 4cr rms ). That technique is often used for high-redshift 
QSOs, and is appropriate if the line originates in the bulge. How¬ 
ever, for disc-dominated line emission such as that seen in this 
source, that approach and the resulting conclusion would be in¬ 
correct. This is a crucial point, given that [CII] emission in other 
high-redshift sources also shows evidence of large-scale ro tation 
dWang et al.ll2013l: IWillott et all2013l: |Pe Breuck et~aI1l2014l) . 

What are the results if we apply the same analysis to other 
QSOs with [CII] detections? Unfortunately, there are almost no 
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sources in the literature with equivalent measurements for perform¬ 
ing this analysis. Most of the QS Os with a [CII] detection do not 
have virial SMBH mass estimates JPe tv et all2 004l;lMaiolino et al.l 


120091: iGallerani et alll2012l : I Wage et aljl20ialWang et al.ll2013h or 

do not show a [CII] velocity gradient (IMaiolino et al. 2005. 120121 : 

IWalter et al.ll2009l : I Venemans et alJl20 12: Cicone et al. 2015j). The 


one source that allows a direct comparison is CFHQS J0210—0456 
Jwillott et alj|2013l) . which has a velocity gradient in the [CII] 
line emission as well as a SMBH mass determined from the 
width of Mg IIA2799; like J1554+1937, its inclination is unknown. 
The SMBH mass is M = 8.0±| g x 10 7 M 0 dWillott et al.l 
120101 ). Using a velocity dispersion of a ~ 80 km s 1 (as es¬ 
timated from the FWHM-u relationship for a Gaussian distribu¬ 
tion), Iwillott et all 120131) found the SMBH in this source to be 
approximately an order of magnitude more massive than predicted 
from the local M-a relation. When we apply our above analysis 
to this source, we obtain the following results. The [CII] line of 
J0210—0456 appears to have a half-width at 20 per cent intensity 
of ~ 140 km s' 1 . Equation [2] suggests a bulge velocity dispersion 
of a ~ 66 km s' 1 /sin(i). If the [CII] line in this source origi¬ 
nates from a rotating disk that has an inclination of 60° to the line 
of sight, we estimate a value for the bulge velocity dispersion that is 
similar to what Willott et al. found by assuming that the line orig¬ 
inates from the bulge itself. For more highly inclined disks, this 
source would be offset by an larger factor than found by Willott et 
al. Thus, when assuming that the [CII] line in J0210—0456 orig¬ 
inates in a rotating disc in a galaxy that als£_contains a bulge, we 
agree with the conclusion o flWillott et al.l d2013l) that the SMBH in 
this source is likely to be over-massive compared to the local M-a 
relation. 

Recently, Iwillott et aid d2015t) updated their discussion of the 
M-a relation for J0210—0456 (and the other known z > 6 QSOs) 
using a similar method to the one presented here, by estimating 
the bulge velocity dispersion from the disc circular rotation veloc¬ 
ity. Their approach and ours are based on the same method, but 
use different formulations; they have applie d the formulation of 
the vo-a correlation presented bv lHol (20071 ). They assign an in¬ 
clination of i = 64° based on the minor/major axis ratio, and find 
a = 98 ± 20 km/s; for the same inclination, our method yields an 
estimate of cr e /2 = 78 km/s. Willott et al. extended their analysis 
to include other z > 6 QSOs in the literature, some of which have 
black hole mass estimates based on the Eddington limit. Where 
inclination estimates were not available, they assumed i — 55°, 
which is appropriate for a population of uniformly/randomly ori¬ 
ented discs. Similar to the results presented here, Willott et al. con¬ 
clude that many z > 6 QSOs may be consistent with the local M-a 
relation, while others have black holes that are much more massive 
than predicted from the local relation. 

Variations in the M-a relation, and particularly its redshift 
evolution, remain controversial (e.g., Kormendy & Bender 2011, 
and references therein). For example, there are significant com¬ 
plexities related to inc onsist encies in velocity dispersion determi- 
nation 1 Bellovary et akOoMh . questionable SMBH ma ss estimates 
(Fine et al.ll2010h . and selection bias at high redshift l iLauer et al.1 
120071) . The analysis of J1554+1937 presented in this paper depends 


on empirical relations between bulge and disc properties, which ne¬ 
cessitates the assumption that this source has a central bulge typical 
of other galaxies. However, the extreme luminosity of this source 
demonstrates that it is atypical, and therefore such relationships 
might not apply. For example, such a luminous source may be ex¬ 
periencing or recovering from a major merger, and therefore would 
not behave like a stable system. 



a [km/s] 


Figure 5. Reproduction of the local M-a relation IY<>m j Kormendy & He] 
1201 3!) (for a = a e / 2 )' Filled circles represent local elliptical galaxies 
and filled triangles represent local spiral galaxies with classical bulges. The 
filled star sy mbols represen t the locations of J1554+-1937 (this paper) and 
J0210—0456 [Willott et alj[2013i) if the value of ry were equal to the ob¬ 
served the velocity dispersion of the [CII] emission line. Hollow star sym¬ 
bols indicate our estimates of velocity dispersion (Sec l4.51 for various in¬ 
clination angles in increments of 15° (as labelled). 


If our assumptions are valid, then we have identified an ex¬ 
tremely luminous, high-redshift QSO that is likely to be consistent 
with the local M-a relation. This result does not contradict the 
conclusions of many earlier observations that high-redshift QSOs 
are generally more massive than the local relation, which have been 
interpreted as indicating that SMBHs in the early Universe grew 
more quickly than their ho st galaxies [Peng e t al. 2006; I Wang et al.l 
1201(1 : Iwillott etall l2013h . Rather, this source provides an addi¬ 
tional data point to probe the extent of this relation’s variation, and 
suggests that some galaxies, like their SMBHs, could have evolved 
quickly at high redshifts. However, before stronger claims can be 
made regarding these conclusions, detailed follow-up of this source 
at higher angular resolution would be valuable, and other extreme 
sources should be similarly analysed. 


5 SUMMARY 

We have presented an ALMA detection of an unusual [CII] line 
in QSO J1554+1937, one of the most luminous QSOs in the ob¬ 
servable Universe, with bolometric luminosity of Lbol = 1.3 x 
10 14 L 0 . The line presented here has a very broad (FWHM ~ 
735 km s' 1 ) non-Gaussian profile, and represents the first [CII] 
detection of this kind. The central velocity of the line indicates a 
redshift of z = 4.627, corresponding to an offset of 800 km s' 1 
from the redshift of broad UV lines, which is typical for QSOs. The 
redshift of the [CII] line seems to be the systemic redshift of the 
system. The velocity map shows a gradient from the east-northeast 
to the west-southwest across the 2 kpc spatial extent of the source. 
However, as the source is barely resolved, it is not currently pos¬ 
sible to determine whether the velocity gradient originates from a 
rotating disc, from bi-directional outflows, or from another origin 
altogether. The profile of the line is consistent with the hypothe- 
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sis of a rotating disc origin, with a maximum circular velocity of 
~ 480 km s' 1 . 

Following the hypothesis that the [CII] line originates from 
a rotating disc, we employed an empirical relation to estimate the 
bulge velocity dispersion in the QSO host galaxy, which requires 
the implicit assumption that the host galaxy contains a classical 
central bulge. The results of this analysis suggest that this source 
is likely consistent with the local M-a relation, or is offset from 
the relation by no more than one order of magnitude in black hole 
mass. This result stands in contrast to earlier results for other high- 
redshift QSO sources, whose black holes tend to be one-to-two or¬ 
ders of magnitude more massive than predicted by this relation.Our 
results are consistent with those of lWillott et al.1 120151) , who used 
a similar approach to update estimates of velocity dispersion for 
z > 6 QSOs. 

Our result stresses the importance of taking into account the 
true dynamical origin of the [CII] line when using it to estimate 
host galaxy parameters. If we were to use its velocity dispersion 
as a proxy for the host galaxy velocity dispersion (equivalent to 
assuming the line originates from a host galaxy bulge), we would 
have reached a very different conclusion about whether this source 
is offset from the local M-a relation. This analysis is important 
because numerous high-redshift sources have [CII] profiles that 
are suggestive of large-scale rotation of the emitting region.s 

A thorough understanding of the [CII] line in J1554+1937 
requires higher-spatial-resolution observations, which are possible 
using an extended configuration of ALMA. ALMA has only re¬ 
cently opened this new window into the high-redshift Universe, 
where high-angular-resolution observations of sub-millimetre lines 
can quickly probe galaxy kinematics at high redshift. Similar ob¬ 
servations of other sources will likely follow, enabling us to deter¬ 
mine whether this type of emission line is truly rare, or if the line 
presented here is simply the first of many to follow. 
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